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Abstract 

Electrodynamics, Continuum Mechanics, Constitutive Theory, Rela- 
tivity Applications of the covariant theory of drive- forms are considered for 
a class of perfectly insulating media. The distinction between the notions 
of 'classical photons' in homogeneous bounded and unbounded stationary 
media and in stationary unbounded magneto-electric media is pointed out 
in the context of the Abraham, Minkowski and symmetrized Minkowski 
electromagnetic stress-energy-momentum tensors. Such notions have led 
to intense debate about the role of these (and other) tensors in describ- 
ing electromagnetic interactions in moving media. In order to address 
some of these issues for material subject to the Minkowski constitutive 
relations, the propagation of harmonic waves through homogeneous and 
inhomogeneous, isotropic plane-faced slabs at rest is first considered. To 
motivate the subsequent analysis on accelerating media two classes of elec- 
tromagnetic modes that solve Maxwell's equations for uniformly rotating 
homogeneous polarizable media are enumerated. Finally it is shown that, 
under the influence of an incident monochromatic, circularly polarized, 
plane electromagnetic wave, the Abraham and symmetrized Minkowski 
tensors induce different time-averaged torques on a uniformly rotating 
materially inhomogeneous dielectric cylinder. We suggest that this ob- 
servation may offer new avenues to explore experimentally the covariant 
electrodynamics of more general accelerating media. 



1 Introduction 

This is paper II of a series of two papers. In paper I it was shown how the 
notion of a Killing vector field on spacetime could be used, together with a 
divergence-less total stress-energy-momentum tensor for a material continuum 
interacting with the electromagnetic field, in order to establish the dynamical 
classical evolution of the medium. It was emphasized that the precise form of an 
electrodynamic force or torque depends on the nature of the decomposition of 
the total stress-energy-momentum into parts describing its electromagnetic in- 
teraction with the medium and the electromagnetic constitutive relations for 
the macroscopic Maxwell equations. The motivation for this approach was 
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to establish whether such quantities could be used to discriminate by exper- 
iment between a number of electromagnetic stress-energy-momentum tensors 
that have been proposed in the past to describe the electromagnetic interaction 
in dielectric media and to provide a comprehensive framework for the analy- 
sis of electrodynamic problems in accelerating media. In this paper, we adopt 
the conventions of paper I and apply the formulation to the computation of a 
time-averaged electromagnetic torque produced in stationary and rotating me- 
dia by an incident plane harmonic electromagnetic wave using the electromag- 
netic stress-energy-momentum tensors proposed by Abraham and Minkowski. 
The first sections demonstrate how the notion of the classical 'photon' in station- 
ary media can be formulated in this framework. It is emphasized that the linear 
momentum of such photons not only depends on the choice of an electromag- 
netic stress-energy-momcntum tensor in the medium but whether the medium 
has an interface with the vacuum or anisotropic and dispersive properties. It is 
then shown how the existence of a family of transverse electromagnetic modes 
satisfying Maxwell's macroscopic equations subject to the Minkowski constitu- 
tive relations for a uniformly rotating homogeneous insulating medium can be 
used to solve the boundary value problem for a plane harmonic electromagnetic 
wave incident on a homogeneous rotating slab. That the time-averaged electro- 
magnetic torque on such a slab is the same for the Abraham and Minkowski 
stress-energy-momentum tensors motivates an analysis of torques on rotating 
inhomogeneous cylinders. In this case we demonstrate that the time-averaged 
torques are significantly different. 

2 Classical 'Photons' in Stationary Homogeneous 
Bounded and Unbounded Isotropic Media 

Although the quantization of the electromagnetic field in a bounded polarizable 
medium is non-trivial the notion of the classical 'photon' has been used to high- 
light the different predictions concerning the linear momentum of light obtained 
by adopting different electromagnetic stress-energy-momentum tensors in such 
a medium. These notions can be defined in terms of time-harmonic classical 
electromagnetic field configurations in homogeneous isotropic stationary media 
and the time-averages of their energy and linear momentum in a fixed spatial 
volume V. 

If a scalar field A(r,t) depends on time t, its time-average over any time 
interval T is 



Furthermore, if A(r, i) is a real spatial p-form, B(r, i) a real spatial q-form in 
any frame U and 




A(r,i) 



Re(.4(r)exp[-iu;t]), B(r,t) 



Re(B(r) exp[-iwt]) , 



where -4(r) and B(r) are complex spatial forms with complex conjugates A(r), B(r), 
then 
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Hence 

( A A B )(r) = ±Re(A(r) AB(r)), (1) 

if T = ^j-. For any bounded spatial 3- form a t in frame U with arbitrary time 
variation, its time-average ( a ) over any finite interval of time T is the 3-form 

(a) = ^J^tdt. (2) 

It follows immediately that if a t is time-periodic, but not necessarily harmonic, 
with period T (i.e. a t — cet+r) then ( a ) = 0. 

If K is a spacelike translational Killing vector field the time-averaged linear 
momentum associated with a drive-form having orthogonal components J^, p u K 
in a volume V relative to U is 

n\V] = - I (P U K )- (3) 
co Jv 

Similarly since U is a timelike translational Killing vector field in Minkowski 
spacetime, the time-averaged energy associated with the same drive-form in a 
volume V is 

£ U kV] = I (P u u)- (4) 
Jv 

This classical energy can be parcelled into N 'energy quanta', each of which 
corresponds to that carried by a harmonic plane- wave quantum with energy fko 
in the vacuum: 

£%\V] = Nhu. (5) 

One then defines the time-averaged if-component of linear momentum associ- 
ated with a classical "photon" in a volume V of the medium to be 

pk ~~^u^>j- (6) 

Different choices of stress-encrgy-momentum tensor for the electromagnetic field 
with the same constitutive relation for the medium will in general give different 
values for p^. In particular for plane harmonic waves propagating in a simple 
unbounded non-accelerating homogeneous, isotropic medium (described by rela- 
tive permittivity e r and relative permeability /i r ) in a direction aligned with K, 
the non-symmetric Minkowski electromagnetic stress-energy-momentum tensor 
yields 

UM Afhuj 

Pk = V) 
co 

where N = y/e r [i r is the refractive index of the medium, while that calculated 
from the Abraham tensor yields 

U,AB /ox 
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Furthermore from the symmetrized Minkowski tensor one finds 



U.SM 1 / U.M , U.AB\ fibJ ( A/" 2 + l\ . . 



Pk 



2c V M 



The underlying origin for these distinctions is that (see tables in Appendix A 
of paper I) although the contributions to the energy density p^ is the same for 
the non-symmetric Minkowski tensor T M , the symmetrized Minkowski tensor 
T SM and the Abraham tensor T AB (since the medium acceleration A = 0) the 
contributions to the corresponding momentum densities p v K are different. 

However even if the medium is stationary, homogeneous and isotropic the 
presence of boundaries can modify these results. Given that all physical me- 
dia do have interfaces with other media including the vacuum the role of these 
interfaces between media with different constitutive properties is relevant in dis- 
cussing physical phenomena. To illustrate this point consider a simple medium 
composed of a stationary plane faced slab of arbitrary finite thickness and con- 
stant relative permittivity e r and relative permeability p, r . If a harmonic circu- 
larly polarized plane wave is incident from the left vacuum half-space normally 
on one plane face it will be partially reflected and transmitted. Suppose it 
propagates through the slab into a half-space composed of a different stationary 
simple homogeneous isotropic medium but with constant relative permittivity 
e" r and relative permeability p r . One can readily compute the electromagnetic 
fields in each region from the electromagnetic junction conditions and hence the 
different classical photon momenta in the three different media. If the slab has 
its plane normals parallel to the direction of propagation K = one finds for 
the different classical photon momenta in the slab: 



U.SM Afprfir (1 + A/" 2 ) U.M ^ 2Af 2 Afp r p, r . . 

Pk = — , p K = — (10) 

C £2^2 + ^2 K c £2_^2 + ^2 

U.AB ^ 2p. r fl r Af 

Pk = — (11) 

K C £2_^2 + ^2^2 

where Af 2 = ? r ^ r . Each of these momenta in the bounded slab depends on the 
properties of the medium outside the slab. Furthermore if the right-hand half 
space is taken to be the vacuum (J\f = l,p r = 1) they are different from the 
momenta above for classical photons in homogeneous unbounded media. 

Analogous results can be obtained by calculating the total time-averaged 
classical angular momentum about some point in V using a spacelike rotational 
Killing vector generating rotations about an arbitrary direction in space. Again 
the momentum densities p v K derived from T M and T AB are different so the 
corresponding classical 'photon' helicities in a finite volume V are different. We 
stress that such 'photon' momenta and 'photon' helicities are strictly classical 
notions calculated from time averages of harmonic plane waves in a stationary 
medium. As such they are not subject to the same conservation laws as genuine 
electromagnetic quanta in the medium. 
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3 Classical 'Photons' in Stationary Unbounded 
Homogeneous Magneto-electric Media 

In dispersive media, constitutive relations between the real spatial fields e u , h u , d u , h' 
are, in general, non-local in spacetime. If the medium is spatially homogenous, 
so that it has no preferred spatial origin, then in Minkowski spacetime it is pos- 
sible to Fourier transform the inertial components of these fields with respect 
to space and time, and work with transformed local constitutive relations. 

For any spatial 1-form a u on spacetime with inertial components a^(r,t), 
define their complex valued Fourier transforms 6q[ (k, uj) by 



where a u {r,t) = a%(r,t) e a , a u {k,oj) = 6^{k,uj)e a , a#(k,w) etandke R 3 . 
Then the source-free macroscopic Maxwell system reduces to 

KAe p (k,u) = u)B u (k,w), KAh^u) = -uD p (k,u) (13) 

where the real propagation wave 1-form K = k • dr. The remaining transformed 
Maxwell equations K A B^k, uj) = and K A ~D u (k, uj) = follow trivially 
from (fig)). It also follows trivially that e u (k,uj) A B y (k,w) = (i.e. e^k,^) is 
perpendicular to b u (k, uj)). Similarly, B^k, w) A K = and D a (k, uj) A K = 0. 

Assume that in any inertial frame U the medium is described in terms of the 
real magneto-electric (1, 1) spatial tensors C de (k, uj), C hb (k, uj), ( db (k, uj), C (k, uj) 
satisfying the symmetry conditions [TJ [5] : 



( d ^(k,uj) = ( de (k,uj), (^(k,uj)=<; hb (k,uj) and ( db \k,uj) = -( he (k,uj), 



where the adjoint of any spacetime tensor T which maps p-forms to p-forms 
is defined by: 

aA*T(/3) = ft A (a) for all spacetime p-forms a, /3 (14) 

and the dispersive magneto-electric constitutive relations in the U frame are 
defined by 



These will (by convolution) give rise to non-local spacetime constitutive rela- 
tions. 



Substituting (|15p and (fT())) in ([TU1) yields a degenerate 1-form eigen-equation 
for e u (k,uj): 

uj 2 ( de (k,uj)(e u (k,uj))+uj( db (k,uj) (#(K A e u (k,uj))) 

+ c# (K A C he (k,w)(e u (k,w))) + # (K A C' ib (k,c) (#(K A e u (k,uj)))) = 0. 

The field b (k, uj) then follows from (fT3"|) . (up to a scaling) and d (k, uj), h (k, uj) 
from ([T5|) . (fT^|) respectively. Equation (fT7|) may be written 




(12) 




(k,oj) 
\k,uj) 



C de (k, uj) (e u (k, uj) ) + C db (k, uj) (b u (k, uj)) 
C he (k, uj)(e u (k, uj)) + C hb (k, uj)(b u (k, uj)). 



(15) 
(16) 



V(k,uj)(e u (k,uj)) = 0, 



(18) 



5 



defining the (1, 1) tensor T>(k, uj). For non-trivial solutions e (k, the deter- 
minant of the matrix T>(k, uj) representing T>(k, uj) must vanish: 

det(2?(k,w)) = 0. (19) 

Note that, in general, the roots of this dispersion relation are not invariant 
under the transformation K — > — K. If one writes k = k|k| in terms of the 
Euclidean norm |k|, and introduces the refractive index JV = |k|^j- > and k 
in place of k, then solutions propagating in the direction described by k with 
angular frequency ui > correspond to roots of (fT9")l (labelled r) that may be 
expressed in the form J\f r = J- r (k,uj). Thus, there can be a set of distinct 
characteristic waves each with its unique refractive index that depends on the 
propagation direction k and frequency uj. When the characteristic equation f| 19[) 
is a quadratic polynomial in Af 2 and has two distinct roots that describe two 
distinct propagating modes for a given uj, the medium is termed birefringent. 
Roots Af 2 such that Af r (k,uj) ^ Af r (— k, uj) imply that harmonic plane waves 
propagating in the opposite directions ±k have different wave speeds. 

Each eigen-wave will have a uniquely defined polarization obtained by solving 
the independent equations in (|18[) for e U ' r (k,uj), up to normalization. Since 
e u,r (k, uj) is complex, it is convenient to introduce the eigen-wave normalization 
by writing 

e u ' r {k,uj) = e u > r (k,uj)n r {k,uj), 

in terms of the complex 0-form e U ' r (k,uj) and complex polarization 1-form 
n r (k, uj), normalised to satisfy 

n'(k,w)A#ir(k,w) = #1 (20) 

for each r. If one applies #e u > r (k, uj) A # to (TT7)) . making use of the symmetries 
between the real magneto-electric tensors C, de (k, uj), (, db (k, uj), C (k, uj), £' lb (k, a;), 
and evaluates it with the eigen-wave e^'^k, uj), one obtains the real 0-form dis- 
persion relation for the characteristic mode r: 

c 2 #(^A#C de (k,c)K)) + ,#(FA#( di (k,,)(#(KAn r ))) + 
+uj#(n F AKA( he {k,uj)(n r )) + # (W A K A C' lb (k, u) (#(K A n r ))) = 0, 

where n r = n r (k, uj) and K = —J\fk ■ dr in terms of Af and k. 

For illustration, consider an unbounded magneto-electric material with 

C de (k,u) = C de (k,w) Id, C' l "(k, w ) = C' l "(k, w ) Id, (21) 

in terms of the scalars ( de (k, uj) X hb (k, ui) and the rank-3 identity tensor Id in 
space. The medium is oriented in the laboratory spatial basis {d x ,d y ,d z } so 
that ( db (k,uj) takes the particular form 

( db (k,uj) = ^i(k,w) dx®d y + /3 2 (k,uj)dy®d x , (22) 

in terms of the real scalars /3i(k, w),/32(k, uj). The matrix representing ( db (k,uj) 
in the laboratory basis takes the form 

/ p 2 (k,oj) \ 
[C db (k,uj)}= ^{k,uj) . (23) 
\ / 
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It follows thalQ 



C he (k,uj) = -/3 2 (k,w) dx®d y - ^(k,uj) dy®8 x . (24) 

Consider a (complexified) harmonic plane wave propagating through the medium 
along the z-axis and polarized in the y-direction: 

e u ' r (k,u) = £exp[ikz - iut] dy (25) 

for some (real) constant £, angular frequency u) > and wave number k. From 
(HHJ), the dispersion relation associated with a polarized eigen-mode e (k, u)) is 

C de (k,uj)uj 2 - C' lb (k,cj)fc 2 - 2/3 2 (k,w)fcw = 0, (26) 

describing propagation in a direction determined by sgn(fc) d z with phase speed 
\oj/k\ depending on the values of C de (k, w), C (k, ui) and /3 2 (k, ui). 

Using (fT5)) and (j2"5")l enables one to construct the 1-forms e u ' r (li, ui), b^'^k, to), 
d^'^k, lu), h U ' r (k,uj). Then from (|6]), the classical photon momenta associated 
with different electromagnetic stress-energy-momentum tensors (see Appendix 
A of paper I) may be calculated. The different magneto-electric photon mo- 
menta associated with the direction K = d z are: 

U,M ( \ fe , U.AB ( \ ^W 2 

Pk,±H = hk ±, Pk.± M = jr 

U.SMt \ 1 ( U,Mf \ i U.AB, A 

Pk,± M = 2 \Pk,±( w ) + Pk,± HJ ' 
where fc± are the roots of 

= c , lfc ^ (~& (k, n) ± ^/ /? 2 2 (k^j + C de (k, c)e b (k, 

with k = (0,0, fc). These momenta are independent of /3i(k, w) owing to the 
choice of wave polarization. For /3 2 (k, w) = one writes C de (k, cj) = e e r (k, w), 
(C' ib (k, w)) _1 = /i Av(k, w) and the magneto-electric photon momenta reduce 
to the classical photon momenta in a stationary, but dispersive, homogeneous, 
isotropic, unbounded polarizable medium, ([7|)-([9]), with a frequency dependent 
refractive index. 



4 Transverse (TDB) Modes in Uniformly Ro- 
tating Homogeneous Media 

In <J3J one used the result that plane harmonic electromagnetic waves can freely 
propagate in a simple homogeneous isotropic non-dispersive non-accelerating 
unbounded polarizable medium. For bounded media one expects that boundary 
conditions will put constraints on the propagation characteristics. For inhomo- 
geneous bounded media single harmonic plane waves are no longer supported 
and if the medium is accelerated (whether bounded or not) finding solutions to 
Maxwell's equations in the medium becomes non-trivial in general. 

iFrom the adjoint relation £ di> t(k, u>) = —( he (k.,uj). 
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A significant difference between the electromagnetic stress-energy-momentum 
tensors advocated by Minkowski and Abraham is that irrespective of the elec- 
tromagnetic constitutive relation describing the medium the Abraham tensor 
depends explicitly on the bulk 4-velocity field of the medium. This in turn 
implies that its divergence will depend explicitly on the bulk 4-acceleration of 
the medium. This characteristic feature is frame independent and gives rise, 
in general, to a non-trivial coupling between the medium acceleration and the 
electromagnetic field. This interaction is of course absent for media at rest or 
moving with constant linear velocity in any inertial frame. By contrast a uni- 
formly rotating medium should be sensitive in principle to such an interaction. 
In the following we will attempt to calculate the significance of this effect. If 
it is possible to measure the torque on a uniformly rotating medium as a func- 
tion of rotation speed this should in principle discriminate between the two 
electromagnetic stress-energy-momentum tensors. 

To expedite this program it proves necessary to excite, by some means, 
electromagnetic fields in an electrically neutral rotating polarizable medium. 
Since in practice one must deal with finite media and fields outside the body 
one is confronted with a difficult problem of electromagnetic scattering from a 
bounded moving polarizable medium. To circumvent this we shall approach the 
problem in terms of the transmission of incident plane harmonic waves through 
a thin cylindrical slab uniformly rotating about its axis of symmetry. If the 
radius of the circular cylinder greatly exceeds its thickness it is reasonable to 
neglect the boundary conditions on its rim. 

Player [3] was one of the first to explore the propagation of waves in a 
uniformly rotating (unbounded) medium using the non-relativistic Minkowski 
constitutive relations. He assumed a rigidly rotating medium with angular speed 
rQ -C Co for all points a distance r from the axis of rotation and concluded that 
in this approximation a certain type of plane harmonic wave in the medium 
was subject to a dispersion relation dependent on f2. Later Gotte, Barnett and 
Padgett [1] argued that there should be a spectrum of such modes. 

We look first for harmonic plane wave modes in a uniformly rotating homoge- 
neous uncharged medium by supposing that the magnetic induction field h u and 
electric displacement field d u are transverse to the propagation of a wave (with 
wave number k and angular frequency u>) along the axis of rotation. In an inertial 
frame U introduce the spatial cylindrical coframe {e 1 = dr,e 2 — rd9,e 3 = d z} 
in cylindrical polar coordinates (r, 0, z) centred at an interface of the medium 
and, for some complex amplitudes B(r), D(r), assume that a pair of the com- 
plexified interior spatial fields take the form 



where m € Z and u> > 0. Such fields will be referred to as circularly polarized 
TDB modes. For constant fl the spatial vector V u = fl-^g describes a rigid rota- 
tion in U about the z-axis in these coordinates and will be adopted to describe 
the bulk motion of the mediunn. Then to first-order (in — ) the constitutive 

v CO ' 

2 With Q constant a real unconstrained medium would not remain unstressed in a strictly 
rigid state according to Newtonian continuum mechanics. Furthermore we eschew all issues 
associated with notions of relativistic rigidity assuming that they are peripheral to the main 
discussion here. 



B(r)exp[i(m + 1)9 + ikz - iwf] (e 1 + ie 2 ) 
£>(r)exp[i(m + 1)9 + ikz - iut] (e 1 + ie 2 ), 



(27) 
(28) 
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relations derived^ from (2.23) of paper I for a simple rotating homogeneous 
medium determine the interior electric and magnetic fields as 



D(r) ( f 

c : exp[i(m + 1)6 + lkz - iwt] [ — — (e 1 + ie 2 ) + 1 

eQ€ r \ e r fi r 

h u = exp[i(m + l)0 + ikz - iut] ( ^^-(e 1 + ie 2 ) - ( f — 



B{r)VLre 3 
D(r)flre 3 



These fields must be compatible with the source free Maxwell system for time 
harmonic fieldfH: 



de u = icuB 1 



dh u = -iLuT) u . 



Substituting into the first Maxwell equation above requires a 2-form to be zero 
(i.e. each component to vanish): 



e'Ae 2 



e'Ae 3 



e 2 A e 3 



dD(r) _ mD(r) 
dr r 

dB{r) B(r) 
dr r 

D(r) = ie e r 



B(r) ikD(r) 

— , tfW 1 

(1 — e r [i r )ll \ r 6q r 



This system has the solution 
D{r) = Ar n 
B(r) = 



uj (to + f )Q 
k k 



ik[i r A 



B(r). (29) 



eo [(m + 1)(1 — e r fj, r )fl + e r /i r Lu] 
for some constant A. Similarly the second Maxwell equation above gives 



e x he 2 



e'Ae 3 



e 2 Ae 3 



with solution 



dB(r) 


mB(r) 


dr 


r 


dD{r) 


D(r) 


dr 


r 


B(r) = 





lo (m + l)fl 
k 



jl r U)- 



D(r) ifc B(r) 



k 



Mo r 
D(r). (30) 



B(r) 
D(r) 



Cr m , 



ike r C 



Ho [(to + 1)(1 — e r /v)fi + e r /j, r ui] 
for some constant C. The solutions (|2"9")l and (fBTJl) are compatible provided 

TV 2 



k TT~ 

Co 



u-(ro+l)fi 1- 



€ r fl r 



(31) 



1 See appendixfor further details of this calculation. 

2 For homogeneous media it is only necessary to solve these coupled equations since the 
other spatial Maxwell equations are then automatically satisfied. 
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For real ui a root k of this dispersion relation may become complex for certain 
values of SI (and describe evanescent waves). Restricting to real k roots, with 
real Q, and real N > such that 



u-nii > o. 

yields propagating solutions with wave numbers as real roots of the dispersion 
relation (f3~Tj). These determine the two possible directions of propagation of 
each TDB mode labelled by m. Let 



kR,L = VR,L 

CO 



(m + 1)^(1 — 



(32) 



where tjr — 1 and t]l = — 1. Then sgn(k^) > and k R denotes a right- 
moving wave while sgn(fcf,) < and kh denotes a left- moving wave. Thus, for 
propagating TDB waves satisfying the dispersion relation (|3ip 

WA+(k R , L ) 

^R,L — —VR,L- 



eo£ r co 



and 



R / ^ iA/"^+(fc fl;L ) m 
BR,L{r) = Vrl r 

D(r) = A+(k R , L )r m 

for the arbitrary constant A^k^^). In this notation the left and right propa- 
gating, left and right circularly polarized, TDB modes in a uniformly rotating 
homogeneous dielectric cylinder can be written: 

e R% = A K%*- L) ex P[i(™ + + ik R,L* - H (e 1 ± ie 2 - ^R,L (l - ^) ^e^) r m 

h U,± = , flii ^^» cx p[ i(m + 1)0 + ifcfl L2 _ iujt ] (e 2 zp ie l )r ™ 

d flL = A^,(k kiL )exp[i(m + l)e + ik RtL z-iojt] (e 1 ±ie 2 )r m 



R,L 

h u,± 

R,L 



A ™ (k ^ L)C0 e x p[i(m + 1)0 + ik RtL z - iut] ( V R,L (e 2 T ie 1 ) - (l - ^) ^e 3 ) r». 



If the axis r — is in the medium and there are no axial sources then the integer 
m > 0. All integers m are permitted if the medium is a shell with the rotation 
axis excluded. Fields with m > are clearly unbounded if r extends indefinitely. 
All modes are eigen- forms of the rotation operator \Cq s with eigenvalue m + 1. 
The interesting physical modes in an unbounded medium are the harmonic plane 
waves with m = 0. These are the waves first explored by Player [3]. One expects 
that these modes could be excited by a harmonic plane wave normally incident 
on a rotating slab from the vacuum. In a simple medium its refractive index Af 
is dispersion-free and independent of the harmonic wave polarization. 

There also exists a set of circular polarized TEH modes that can be calculated 
using an initial ansatz where the magnetic field h u and electric field e u are 
transverse to the direction of propagation along the rotation axis: 



e 



E(r)exp[i(m + 1)6 + ikz - iwt] (e 1 ± ie 2 ) (33) 
#(r)exp[i(m + 1)0 + ikz - iuit] (e 1 ± ie 2 ) (34) 
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with to s Z as before. However if 7^ it appears that there are no harmonic 
plane wave TEH modes among this set, even for m = 0, so we will not discuss 
them further here. 

If a plane circularly polarized harmonic wave in the vacuum is incident nor- 
mally on a homogeneous uniformly rotating thin cylindrical slab consisting of 
material satisfying a simple constitutive relation (with e r , [i r constant) then it is 
straightforward to satisfy the electromagnetic jump conditions at the slab plane 
interfaces with left and right moving m — 0, circularly polarized, TDB modes 
in the medium (and a reflected and transmitted circularly polarized, harmonic 
plane wave in the vacuum). From equation (3.2) of paper I, one can use such 
a solution to calculate the real instantaneous electromagnetic K-drive on the 
slab. Such a drive will fluctuate harmonically in time with angular frequency 
a multiple of the incident frequency. For any t^ em the field in the slab is 
time- harmonic in the laboratory frame U and the time-average of C\j p^ EM is 
zero. Furthermore since J K ' E is a real quadratic function of the time-harmonic 
fields e u ,h u ,d u ,h u in the medium, the time-average of JjJ E is independent 
of the propagation coordinate z. Hence from equation (3.2) of paper I, the 
time-averaged total i^-drive on a fixed volume V of the rotating medium is: 



The surface dV consists of the two opposite interior circular discs and the interior 
cylindrical rim. For all the tensors in Appendix A of paper I, the contribution to 
the integral from the latter is zero since the integrand in ([33)) is zero when pulled 
back to the interior rim. Furthermore the contributions to the integral from the 
two opposite interior faces cancel since the time-averaged ( j^ m - EM ) j s inde- 
pendent of z. Hence one concludes that if a normally incident plane harmonic 
vacuum wave excites a plane harmonic TDB mode in such a simple (uniformly 
rotating or stationary) medium the induced time-averaged if-drive will be zero. 
Thus one cannot discriminate between Minkowski and Abraham stress-energy- 
momentum tensors by exciting plane TDB modes in a simple medium in this 
manner. One alternative is to excite harmonic fields in a medium in which the 
time-averaged ( J^' in,E ^ j s no t independent of z. Then the contribution from 
opposite interior faces will not necessarily vanish. The simplest strategy to im- 
plement this requirement is to employ fields in an inhomogeneous dispersion-free 
medium in which the permittivity or permeability varies with z along the axis of 
the slab. Unfortunately the simple TDB plane wave modes above are no longer 
solutions to the Maxwell equations in such an inhomogeneous medium. Never- 
theless it will be shown below how the computation of the Maxwell fields in such 
cases can lead to a method of discriminating between different electromagnetic 
stress-energy-momentum tensors in uniformly rotating inhomogeneous media. 

5 Plane Waves Incident on a Stationary Inho- 
mogeneous Dielectric Slab 

The theory of Killing forms was exploited above in order to calculate the time- 
averaged harmonic electromagnetic forces and torques on perfectly insulating 




11 



homogeneous polarizable media described by the Minkowski constitutive rela- 
tions. For stationary or non-relativistically uniformly rotating media it was 
found that the induced electromagnetic drives did not discriminate between the 
electromagnetic stress-energy-momentum tensors proposed by Minkowski and 
Abraham. If one stays with slabs excited by incident plane harmonic vacuum 
electromagnetic waves it is natural to enquire about the influence of material 
inhomogeneities on this result. However a polarized plane harmonic wave inci- 
dent in the vacuum on a stationary inhomogeneous slab of polarizable material 
will not in general propagate in the medium as a polarized plane harmonic wave. 
To facilitate a discussion in the next section of the behaviour of such an inci- 
dent wave on a rotating slab, a solution to the problem of fields in a stationary 
inhomogeneous slab will be considered in this section. 

The stationary slab is oriented in a Cartesian frame with coordinates {x, y, z} 
as indicated in figure [TJ Suppose a plane polarized plane wave monochromatic 
wave with angular frequency u) > 0, propagating in the z-direction, is incident 
on a slab of thickness I and infinite extent in the x and y directions. 

Let the slab have relative permittivity e r (z) and constant relative permeabil- 
ity [i r . The slab has parallel interfaces (with the vacuum) at z = and z = i 
thereby distinguishing the spatial region z < denoted I, < z < £ denoted 77 
and z > I denoted III. 

x 

Incident Wire 



KM mode Transmit ted Wave 



Bfiflected Wive 



/ n in 



j = z=t 



Figure 1: Geometry of the inhomogeneous stationary dielectric slab and the 
propagation direction of waves in the vacuum regions. 

Let F 11 = dA 11 in region 77 with 

A 11 = A A(z)exp[-iut] dy 
for some (complex) function A{z) and (real) constant Ao with physical dimen- 
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sions (so that A(z) is dimensionless) . Then 

F u = dA u = dz Ady- iuA(z)dt A dyj A exp[-iut] . 

The medium is at rest relative to the inertial frame U so is assigned the 4- 
velocity V — —dt- Thus, from the covariant dielectric constitutive relation 
(2.23) of paper % 

, n _ ( 1 dA(z) 



G = I — dz A dy — ie r (z)uiA(z)dt A dy I Aq e o exp[— \ujt] . 

For an uncharged medium (J = 0) the inhomogeneous Maxwell equation (2.21) 
of paper I yields 

, ^rr ( 1 d 2 A(z) u 2 e r (z)A(z)\ A r . . , 

d*G u = H , Mo coeo exp -mf dzAdxAdt = 0. 

If the relative permittivity is taken to be the function, e r (z) = a + 4r for some 
(dimensionless) constants a, /?, this Schrddinger-like differential equation can be 
solved in terms of a basis of Airy functions of the first and second kind (Ai and 
Bi respectively) to give 

A{z) = CiAi[/C(*)] + C a Bi[/C(z)] , 
where Ci , C*2 G C and 

i 

(at \ ( Lu 2 j3^ r 



f3 , ~J\ lc 2 

Thus, for this inhomogeneous dielectric with linearly-varying permittivity: 

A 11 = A ( CiAi[/C(z)] + C 2 Bi[lC(z)] S J exp[-iwt] dy. (36) 

The dimensionless complex constants C\ , C*2 can be determined from two inter- 
face conditions. Consider vacuum electromagnetic plane waves in region / with 
(real) amplitude Aq and complex amplitude E^, propagating in the (positive 
and negative) z-directions respectively: 

A 1 = Aq ( exp[ifc^z — iiot] + _E£exp[ifc£z — kjt\\ dy, (37) 

where kn and k^ denote distinct roots of the vacuum dispersion relation 

k 2 -^- = (38) 
c o 

with sgn(fc|j) > and sgn(fc£) < 0. Similarly, a transmitted field in region /// 
is written 

A 111 = AoE^'e^ik^z-iuj^dy (39) 

for some complex amplitude E^ 1 and sgn(fc^ 7/ ) > 0. At each interface one 
must satisfy the interface conditions on a spacetime hypersurfac^f). On the left 



3 See appcndixfor their formulation on spacetime. 
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interface / = z and on the right interface / = z — i. If fig (fit) denote the pull- 
back of forms to z = (z = £) respectively, the interface boundary conditions 
become 



n* (F 1 - F 



n: (f h - F 111 ) 







CIq (-kG 1 - *G TI ) = ST,(*G 



II 



G ) = 0. 



These yield a linear system of equations for the dimensionless complex variables 
El,Ci,C2, E^ 1 : 

( „ „ V\uC\ 



M[K(t)] Bi[K(£)] 

Ai'[/C(0] Bi'[/C(£)] 

uv— Ai'[/C(0)] Bi'[/C(0)] 

V -1 Ai[/C(0)] Bi[/C(0)] 



-exp 



Co 



-i/i r — exp 
co 






co 



( E L \ ( 

Ci 



1//,. 



V i 



co 



where 



Ai' [AC(c 



Ai[/C(z)] 



and similarly for Bi'[«;(a)]. The system admits a solution provided the 4x4 
complex matrix above is non-singular. 




Figure 2: Electric field amplitude demonstrating transmission of a normally 
incident harmonic plane wave from the left through a stationary linearly- 
inhomogeneous plane dielectric slab. 



The solution is readily obtained but its complicated structure will not be 
displayed here. It replaces the plane wave dispersion relation in a homogeneous 
medium. Suffice to add that with the aid of (|5d]) - (I5^1) . one may now explicitly 
construct the Maxwell and excitation 2-forms in all three regions, and hence 
the electromagnetic field 1-forms {e u ,h u ,d u ,h u } in all three regions. The 
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schematic behaviour of the real electric electric field amplitude in all regions 
is shown in figure [5] Solutions with different incident polarizations follow in a 
similar manner. 

Once one has a complete Maxwell solution for a given incident wave this 
can be substituted into any if-drive in order to calculate time-averaged force 
and torque pressures on any area of the slab in terms of the incident wave 
parameters Aq,lo and the medium characteristics e r (z),£. For the linear po- 
larized plane waves (|37[) and (|39l) one finds non-zero time-averaged integrated 
electromagnetic torques in any finite volume of the medium associated with the 
Minkowski, symmetrized Minkowski and Abraham drive forms. However they 
are all the same. This is to be expected since, if one observes from the table 
in Appendix A of paper I that with the medium at rest (U — V), although the 
medium is inhomogeneous, the associated current 2-forms for U (K) = arc 
the same and ( Px ) = f° r harmonic waves. Moreover if one considers inci- 
dent circularly polarized harmonic plane waves the total time-averaged torques 
are each zero. Hence harmonic plane waves cannot distinguish the effects of 
the Minkowski, symmetrized Minkowski and Abraham electromagnetic stress- 
energy-momentum tensors when incident on non- accelerating inhomogeneous 
slabs of simple media. It is natural to enquire how this result may change when 
the medium accelerates. This is explored in the next section for a simple planar 
inhomogeneous dielectric slab in uniform rotation about its axis of symmetry. 



6 Circularly-Polarized Plane Waves Incident on 
a Simple Rotating Inhomogeneous Dielectric 
Cylinder 

In the last section the fields in a stationary linearly inhomogeneous slab ex- 
cited by an incident plane harmonic wave were found. Such fields yield a non- 
zero time-averaged net force (and torque) on such a slab that is the same for 
both the Abraham and symmetrized Minkowski electromagnetic stress-energy- 
momentum tensors. One could attempt to find the fields excited by an incident 
electromagnetic pulse such as that produced by a laser. Since such fields are not 
harmonic they may yield forces and torques that distinguish between such ten- 
sors particularly if one also employs a slab with anisotropic or magneto-electric 
properties [5] . However the precise nature of the fields in a laser pulse with a 
finite spot size is clearly more difficult to ascertain theoretically and this makes 
the problem of matching fields at the slab interfaces more difficult than with 
harmonic plane waves. 

In this section it is shown how circularly polarized plane harmonic waves 
can excite a transverse plane wave field in a simple inhomogeneous uniformly 
rotating medium analogous to that found in the previous section. The form of 
this wave is determined by a single amplitude function of z in the cylindrical 
coordinates used in 21 satisfying a second order ordinary differential equation 
with coefficients dependent on the frequency of the incident wave, the speed of 
rotation of the slab and its constitutive properties. Although this equation has 
no solution in terms of simple analytic functions it is amenable to numerical 
analysis. Furthermore one can match a basis of such solutions to the vacuum 
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plane waves at the rotating plane slab interfaces. If one accepts that the edge 
effects produced by the rim of the rotating slab are ignorable then an approx- 
imate numerical estimate of the fields excited in the slab by normally incident 
harmonic plane waves can be made. From such a solution one may compute the 
net time-averaged torque on the slab for both the Abraham and symmetrized 
Minkowski electromagnetic stress-energy-momentum tensors. It will be shown 
that for a simple medium with relative permeability one and a permittivity that 
varies linearly with z in the slab these two torques are significantly different. 

In cylindrical coordinates, suppose the cylindrical slab rotates with constant 
angular speed f2 about its axis of symmetry along the z-direction. Let the 
cylinder domain II have constant permeability /i,. but relative permittivity e r (z) 
and denote the left and right vacuum domains by I and /// respectively. Since 
the field in this medium is to be excited by a transverse circularly polarized 
plane harmonic vacuum wave we look for monochromatic time-harmonic interior 
solutions with a TDB-type field ansatz in the cylindrical coframe used in SjU 

= ccT 1 £B(z)exp[iO - icjt] (e 1 + ie 2 ) 
= e £ D(z)exp[i9 - iuit] (e 1 + ie 2 ). 

With the real constant £ assigned the physical dimensions of rjqj^ , the complex 
amplitudes B(z),D(z) to be determined are dimensionless. The electric field 
1-form eYj and the magnetic field 1-form are now determined from the 
constitutive relation for a rotating dielectric. Generalizing the m — 0, TDB 
mode in homogeneous media found in SjH they may be written in the form. 

efj = £ (E 1 (r,z)e 1 + E 2 (r, z)e 2 + E 3 (r,z)e 3 ) exp[i6» - iw <] 

hYj = c e £ (#i(r, z)e 1 + H 2 (r, z)e 2 + H 3 (r, z)e 3 ) exp[i6» - iut] . 

Using the constitutive relations (2.23) of paper I with V u — Vide it follows that: 

, ^ D(z) , , B(z) 

e r {z) (i r 

E 2 {r,z)=iE 1 {r,z) H 2 (r, z) = Lffi(r, z) 

rB{z){l-e r (z)fi r )n r_D{z) (e r (z)^ r - 1) Q 

c a fi r e r (z) Cf)fj, r e r (z) 

to leading order in ^ . The fields F 11 and G 11 can now be constructed to satisfy 
dF 11 — and d -kG 11 = 0. These equations require 

D(z)Y B(z) 



e r (z) J co(j, r e r (z) 
D(z) 



^ e r (z)fj, r (fl — lj) — Clj 

B'(z) - — ^L(e r (z)^ r (Q-ij)-n), 
c e r (z) V / 

where f'(z) denotes the derivative of f(z) with respect to z. Substituting D{z) 
into the first differential equation gives a second order differential equation for 
B(z): 

B _ yL r (n-wy r (z)B'(z) + (e r {z)y r (Vl-u)-ti) 2 B{z) = q 
e r (z)fi r (Cl — uj) — il CQ[J, r E r (z) 
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Note, that with a real e r (z), the coefficients of this differential equation are real. 
Write the solution for B(z) as 



B(z) 



(41) 



in terms of complex dimensionless constants C\, C 2 and a basis of real solutions 
{(Ti(z), (72(2)}. It follows that 



D(z) 



cpe r (z) (da'^z) + C 2 a' 2 (z)) 
e r (z)fi r (fl — ui) — SI 



The forms F 11 and G 7/ inside the rotating slab are now expressed in terms of 
a basis of solutions for (14TJ1) as: 



F 



11 



co(Cia[(z) + C 2 a' 2 (z)) geMm &0 ^ + 

e r (z)/i r (S2 — ui) — ll ' 
+ (Cicri(z) + C 2 cj 2 {z)) £exp[i0 - Suit] (ie 1 - e 2 ) A e 3 



rf2 (1 - e r (z)^ r ) (Ci<ti(z) + C 2 a 2 (z)) 



G 



11 



c e r {z)[i r 
e e r (z)c (CiCTi(z) + C 2 ct' 2 (z)) 



£ exp[i0 - iwi] e° A e 3 



(42) 



e r (z)/i r (51 — a;) — SI 



£exp[i6» - iwt] e° A (e 1 + ie 2 ) 



£0 
//.,. 



+ — (Ci 0-1(0) + C 2 cr 2 (z))£cxp[i<9 - iu;t] (ie 1 - e 2 ) A e 3 



eprfl (1 - e r (z)n r ) (C^jz) + C 2 a' 2 {z)) 
fj, r ( e r (z)fi r (fl — ui) — DA 



£exp[i# — Suit] e 1 A e 2 



As we have seen in the last section, when SI = and e r (z) varies linearly with 
z, the basis solutions can be expressed in terms of Airy functions. The complex 
constants in (|4"2"j) are determined as before by matching F 11 and G 11 to harmonic 
plane waves at the interfaces of region II with regions I and III. Thus with 
F 1 = dA 1 , F IH = dA 111 , G 1 = eqF 1 , G hi = e F IH and 



A 1 = 



A 



in 



c £ 



c Q £ 



exp[i6* — hit] (exp[ikftz] + i?£exp[ifc£z]) (e 1 + ie 2 
J B^ r/ exp[i(9 - Suit + Sk^z] (e 1 + ie 2 ) , 



the plane interface conditions yield the following equations for the dimensionless 
complex variables E I L ,C\,C 2l E^ 1 : 



I ! oj(0) a 2 (0) 







V 

where 



A(0) 
<n(0) 

X(l) 



A(0) 

(7 2 (0) 

gjCO 

AW 



A(z) 







exp 



iu^ 



CO 



-i^ r exp 



iw^ 
Co 



C 2 

V ^ / 



/ - 1 \ 

lfi r 


V J 



e r (z)jjL r (S7 — ui) — SI 



CO 
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This matrix equation admits a solution provided 



exp 



1 — i 

co 



(( 



, , ^(ai(£)4(0)-a( (0)^(1) 

CTi(£)<T 2 (0) - CTl(0)<T 2 («) 



A(0)A(£) 



. (^ffl^W-^iW^ffl) . (^(OMl)- tnffl^O)) ^ 

aw 1Mr X(o) J 

Thus one has a matched solution to the complete electromagnetic system in 
terms of a basis of solutions for (|40[) and e r (z). With K = de, one may now 
express the time-averaged total torques, — J v (iu dr^ B ) and — f v {iu dr^ M ), 
on a thin rotating disk of volume V in terms of this solution and for each Q and 
uj numerically integrate (|40l) to find the needed values of o~\{z) and oi{z) in the 
slab. 

The results in figures [3] and 2] are for e r (z) = a + j z with real constants a, /3 
describing the variation of permittivity across the cylindrical disc of thickness 
I. They clearly indicate that if one could devise an experiment to measure 
the mechanical torque needed to maintain a simple inhomogeneous dielectric 
medium in uniform rotation both in the absence and presence of a normally 
incident circularly polarized harmonic plane wave then one may be able to 
discriminate between the effects induced by different electromagnetic stress- 
energy-momentum tensors. 



7 Conclusions 

Unlike the Minkowski tensor and its symmetrized version, the divergence of the 
Abraham tensor depends explicitly on the bulk acceleration field of the medium. 
Given simple electromagnetic constitutive relations this implies the existence of 
local acceleration dependent Abraham forces and torques in the presence of 
electromagnetic fields. One could write the Abraham tensor as the sum of any 
other symmetric tensor (including the symmetrized Minkowski tensor) and the 
difference between the two. This of course makes no difference to the divergence 
but assigning the difference in the tensors to the accompanying matter tensor 
appears to us ad hoc and unwarranted. As an illustration we would prefer to 
identify the radiation reaction force experienced by a radiating point charge 
(that depends on the rate of change of the acceleration of the particle) with 
an intrinsic property of the charge's interaction with the electromagnetic field 
rather than with the charge's inertia (that conventionally depends only on the 
acceleration) and some additional force. Furthermore from an action variational 
viewpoint it is natural to associate particular stress-energy-momentum tensors 
with particular interacting systems since the extrema of the total action de- 
termine both the total stress-energy-momentum and the field equations for the 
fully coupled system. 

The view advocated here is that it does make sense to try and pin down by 
theoretical modelling and experiment the most appropriate form of an electro- 
magnetic stress-energy-momentum tensor for as large a range of electromagnetic 
constitutive relations as possible. As emphasized in [2d], we believe that vari- 
ational principles including gravitation are a valuable theoretical guide in this 
endeavor. On the experimental side it is suggested here that experiments in- 
volving accelerating media might offer a sensitive means to discriminate between 
alternatives. 



18 



From the calculations for the macroscopic time-averaged specific torques on 
simple media shown in figures [3] and HI it is clear that those corresponding to the 
symmetrized Minkowski tensor are much larger than those associated with the 
Abraham tensor. In each figure the physical torque has been calculated for a 
cylinder of radius 1 m and thickness I — 1 mm as a function of the uniform rota- 
tion speed ft in radians per second. The parameters defining the linear variation 
of permittivity across the thickness of the cylinder are a = 1 and f3 = 100. The 
torque curves are calculated for a range of frequencies (indicated in Hz in the 
legend on the right of each figure) of the normally incident circularly polarized 
incident plane wave. The magnitude of each physical torque (in Newton m) is 
obtained by multiplying the indicated specific torques by £ 2 where £ is the am- 
plitude of the electric field of the incident plane wave in volts m -1 . Recall that 
the time-averaged power in a harmonic plane wave crossing unit normal area 
is £ 2 /(2Zq) watts m~~ 2 in terms of the characteristic impedance Zq = ywAo 
of free space. The detection of such a torque would require the measurement 
of an externally applied torque to maintain a non-deformed cylinder in uniform 
rotation. It is interesting to note from the figures the torque asymmetry about 
n = Q. This asymmetry can be shifted along the horizontal axes by using an 
oppositely polarized incident beam and is clearly a direct manifestation of its in- 
trinsic angular momentum. The significant quantitative differences between the 
time-averaged Abraham and symmetrized Minkowski torques in simple inhomo- 
geneous rotating dielectrics suggest that new opportunities exist for exploring 
the effects of classical electrodynamics in more general rotating media. 
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A Notation 

The natural mathematical language to discuss the differential properties of ten- 
sor fields on spacetime and their relation to integrals over material domains is 
in terms of differential forms and their associated exterior calculus (Benn & 
Tucker 1988). In this section a brief summary is given of the relevant notation 
used in subsequent sections. A key concept throughout involves the role of the 
spacetime metric tensor field and possible isometries that it may possess. A 
spacetime metric tensor field g is a symmetric bilinear form on spacetime that 
can always be represented in a local cobasis of differential 1-forms {e a } as 

g = -e° <g> e° + e 1 ® e 1 + e 2 ® e 2 + e 3 <g> e 3 . (43) 

If {Xt} is the dual local basis of vector fields on spacetime defined so that 
e a (Xt) = 5£ (a, 6, = 0, 1, 2, 3) one has the induced contravariant metric 

g^ 1 = -A <g> X Q + Xx® Xi+ X 2 ® X 2 + A 3 (g) X 3 . (44) 

If a is any given 1-form, one has an associated vector field a defined so that 
(3(a) = g~ 1 (a,/3) for all 1-forms /3. Since g^ 1 is symmetric this will be abbre- 
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Angular speed (radjsee) 



Figure 3: Variation of time-averaged specific torques with angular speed of a 
thin rotating cylinder calculated with the Abraham electromagnetic stress- 
energy-momentum tensor. The physical torque is obtained by multiplying each 
specific torque by £ 2 where £ is the electric field amplitude in volts mT 1 of 
the normally incident plane circularly polarized harmonic wave with frequency 
(including the optical range) given by the legend on the right of the figure. 




Figure 4: Variation of time-averaged specific torques with angular speed of a 
thin rotating cylinder calculated with the symmetrized Minkowski electro- 
magnetic stress-energy-momentum tensor. The physical torque is obtained by 
multiplying each specific torque by £ 2 where £ is the electric field amplitude 
in volts mT 1 of the normally incident plane circularly polarized harmonic wave 
with frequency (including the optical range) given by the legend on the right of 
the figure. 
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viated a = g 1 {a). In a similar way if X is any given vector field one has an 
associated 1-form X = g(X). 

In the g-orthonormal basis {e a } one has a canonical local 4-form denoted *1 
and defined to be e° A e 1 A e 2 A e 3 . The Hodge map * induced by *1 maps p-forms 
to (4— p)-forms on spacetime (Benn & Tucker 1988; see also Appendix A of paper 
I) . The metric also uniquely defines the covariant derivative Vx with respect to 
any vector field X. This has the property Vx.9 = and \7 X Y — \7 Y X = [X, Y] 
for all vector fields X,Y. In this expression [X, Y] denotes the commutator 
bracket. While Vx has a type-preserving action on any tensor field the exterior 
derivative d is defined to act only on antisymmetric tensor fields (differential 
forms) and has the property do d = 0. Contraction of any p-form j3 with X is 
denoted ixP- The interior operator ix is a graded derivation defined so that 

ix(aAP) = (i x a) A0+(-l) p aAi x /3, (45) 

for any p-form a and g-form j3. If p = 1, one defines ixa — a(X), and if p = 0, 
ixce = 0, for all vector fields X. One has the useful relations 

**$ = (-l) p+1 $ 
i x *& = *($>AX) 
Cx$ = ixd$ + dix&, 

for any p-form $ on spacetime where Cx denotes Lie differentiation with respect 
to X. 

The spacetime divergence operator V- takes a simple form if one uses the 
(/-orthonormal basis above. Acting on a symmetric covariant tensor T it defines 
by contraction on the first argument the 1-form 

3 

V-T = ^(V Xa T)(A a ,-) (46) 

where {A a } = {— Xq, X\, X2, X3}. For a symmetric tensor T it is sufficient 
to write the right hand side as (Vx a T)(X a ). A local spacetime isometry with 
respect to g is a local diffcomorphism that preserves this metric. A vector field 
K that generates such a diffcomorphism is called a Killing vector field and it 
satisfies Cxg = in terms of the operation of Lie differentiation with respect 
to K. The operators *, d, ix, Vx, &x offer a powerful computational tool-kit 
when working with differential forms. 

For any smooth p-form <!> in a bounded regular region M. of a manifold one 
can express the integral of ci$ over M in terms of the integral of $ over the 
boundary dM of A4: 

[ = f (47) 

J M JdM 

This is a statement of Stokes' theorem for p- forms. 

The Gibbs calculus of vector fields in 3-dimensional Euclidean space is read- 
ily exposed by correspondences induced by the exterior operations above. A 
space of 3-dimensions may be considered as a particular hypersurface in space- 
time. If the metric above induces the metric 

g = e 1 O e 1 + e 2 ® e 2 + e 3 ® e 3 (48) 
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on this hypersurface it is Euclidean and one may introduce the Euclidean canon- 
ical form #1 = e 1 A e 2 A e 3 by restriction. Since spacetime is assumed time- 
oriented one may employ a future-pointing timelike unit vector field U on space- 
time (with g(U, U) = — 1 and U = e°) to fix a coherent orientation by relating 
#1 to *1 by 

★ 1 = £/A#l. (49) 

If the hypersurface is given a.s t — constant for some time coordinate t, an 
inertial frame exists in Minkowski spacetime with U = ^ Jj = — Xq such that 
V U = 0. Throughout this article the constant Co denotes the speed of light in 
the vacuum. Any p-form a on spacetime is termed spatial with respect to such 
a U if ijja = 0. An over-dot will denote (Lie) differentiation with respect to the 
coordinate t so 

a = caCjja — Cg t a. (50) 

On a Euclidean hypersurface in spacetime, exterior differentiation of spatial 
p-forms (j> is denoted d<j) such that 

d<p = d<j) - U A £u(f> = d(j) + dt A Ca t <t>. (51) 

One has the following relations 

#1 = -% v *1 = -*U 
##</> = 4> 

for all spatial p-forms 4>. Furthermore if v, w denote Euclidean vector fields 
in the Gibbs notation corresponding to the vector fields v, w for some spatial 
1-forms v,w then 

div v corresponds to #<L# V > (52) 

curl v corresponds to ; (53) 

v x w corresponds to #(u A w) ; (54) 

grad ip corresponds to dtp, (55) 

for any scalar ij) field on spacetime. 



B Electromagnetic Fields on Spacetime 

We suppose that a polarizable material continuum is given in terms of a set of 
piecewise smooth material properties that determine its interaction with classi- 
cal gravitational and electromagnetic fields. The classical macroscopic Maxwell 
system for the electromagnetic 2-form F on spacetime can be written as [5]: 

dF = and d*G=j, (56) 
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where the excitation 2-form G depends on the interaction with the medium and 
the 3-form electric 4-current j encodes the electric charge and current sourc^. 
Such an electric 4-current describes both (mobile) electric charge and effective 
(Ohmic) currents in a conducting medium. To close this system in a background 
gravitational field, electromagnetic constitutive relations relating G and j to F 
are necessary. 

The history of a particular observer field in spacetime is associated with an 
arbitrary unit future-pointing timelike 4-velocity vector field U. The field U 
may be used to describe an observer frame on spacetime and its integral curves 
model idealized observers. An orthogonal decomposition of F with respect to 
any observer field U gives rise to a pair of spatial 1-forms on spacetime. The 
1-form spatial electric field e 11 and 1-form spatial magnetic induction field h u 
associated with F are defined with respect to an observer field U by 

e u = i v F and c b u = i v * F. (57) 

Since g(U, U) = -1 and i v e u = i v h u = 0: 

F = e u AU -*(cob u AU). (58) 

Likewise the 1-form spatial displacement field d u and the 1-form spatial mag- 
netic field h u associated with G are defined with respect to U by 

h u 

d u = i v G and — =iu*G, (59) 
Co 



,'h u 



so 

G = d u A U - * A Uj , (60) 

with iud u — ijjii u = 0. At the history of any sharp interface between different 
media, given as the piecewise smooth (non-null) spacetime hypersurface / = 0, 
the system of Maxwell equations is supplemented by interface conditions on the 
fields F and G 



[F] 
[*G\ 



Adf = 

f=o 

Ad/ =j s , 

f=Q 



(61) 



where [H] denotes the discontinuity in the field H across the hypersurface. The 
3-form j s on the hypersurface is non zero if it supports a real current 3-form 
there. 



4 All electromagnetic tensors in this article have dimensions constructed from the SI di- 
mensions [M], [L], [T], [Q] where [Q] has the unit of the Coulomb in this system. We adopt 
[g] = [L 2 ] , [G] = [j] = [Q] , [F] = M where the permittivity of free space eo h a s the dimen- 
sions [Q 2 T 2 M~ 1 L~3] and cq = ^-1— _ denotes the speed of light in vacuo. Note that the 

operators d and V preserve the physical dimensions of tensor fields but with [g] = [L 2 ], for 
p-forms a in 4 dimensions, one has [*a] = [a][L 4— 2p ]. 
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C Time-Dependent Maxwell Systems in 3-Space 



The spacetime description above is natural for the Maxwell system since it makes 
no reference to any particular frame in spacetime. However to make contact with 
descriptions in particular frames or non-relativistic formulations a reduction in 
terms of frame dependent fields becomes mandatory. The spacetime Maxwell 
system can now be reduced to a family of parameterized exterior systems on M 3 . 
Each member is an exterior system involving forms on R 3 depending paramet- 
rically on some time coordinate t associated with U. Let the (3 + 1) split of the 
electric 4-current 3-form with respect to a foliation of spacetime by spacelike 
hypersurfaces with constant t be 

7 U ~ 

3 = —AU + p u , (62) 
Co 

with iuJ u = iuP U — and p u = p u #l, where J u ,p u are the spatial electric 
current density 2-form and spatial electric charge density 3-form respectively. 
The differential operator d on spacetime forms is adapted to those spacetimes 
(such as Minkowski spacetime where gravity is absent) that can be foliated by 
hypersurfaces with constant coordinate t where U = ^c\. Then, from (|56p . 

dj = (63) 



yields 



dj u +p u = 0. (64) 



The (3 + 1) split of the spacetime covariant Maxwell equations (151))) with respect 
to U = — codt becomes 

de u = -B u (65) 
dB u = (66) 
dh u = J u + t> u (67) 
dT> u = p u (68) 

where T> u = #d c/ and H u = #b c/ . All p-forms (p > 0) in these equations are 
independent of dt, but have components that may depend parametrically on t. 



D Electromagnetic Constitutive Relations 

The two 2-forms F and G in the macroscopic Maxwell equations on spacetime 
are fundamentally related by smoothing the microscopic sources of the elec- 
tromagnetic fields in the medium. In many circumstances one then relies on 
phenomenological relations for closure relations. In such relations the excita- 
tion form G is in general a functional (possibly non-local in spacetime) of the 
Maxwell form F, its covariant differentials, thermodynamic properties, defor- 
mation properties and the state of motion of the medium: 

G = Z[F,VF...]. 
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Such a functional may induce non-linear and non-local relations between d u , h u 
and e u , h u . Electrostriction and magnetostriction arise from the dependence of 
Z on the deformation tensor of the medium and its covariant derivatives. For 
general linear continua, a knowledge of a collection of constitutive tensor fields 
on spacetime may suffice so that 

N 

G = ^^'[VF,...]. 

r=0 

In idealized (non-dispersive) simple continua, one adopts the idealized local re- 
lation 

G = Z(F), 

for some degree 4 constitutive tensor field Z, parameterized by scalars that 
depend on the medium. In the vacuum G — eoF where eo is the constant 
permittivity of the vacuum. Regular lossless, non-conducting, linear isotropic 
media can be described by a bulk 4-velocity field V of the medium, a real 
relative permittivity scalar field e r > and a real relative permeability scalar 
field fi r > 0. In this case, the structure of the tensor Z follows from 

G = e Q e r i v FAV- — * (i v *FAV) 
= e (e r ) i v F A V + —F. 

V Mr / Mr 

For inhomogeneous media the relative permittivity and permeability scalars e r 
and \i r will not be constants. In a general frame U comoving with the medium 
{U = V), (69| yields 

d v = e e r e v and h v = (^Mr)" 1 ^, (70) 

which are the familiar closure relations for simple (idealized) electrically neutral 
isotropic non-dispersive polarizable media. 



D.l Minkowski Constitutive Relations for Moving Media 



The above algebraic constitutive relation (|69|) involves the bulk 4-velocity field 
of a simple medium. It is straightforward to find the induced relations between 
the fields {e u , h u ,d u , h u } relative to an observer in a frame U V. From ([55)1 
and (pDl . one has 

e v = i v F = e u {V)U ~U{V)e u -*(c h u A A ^) 
c h v = i v *F ^★(e l/ AU AV) + c b u (V)U - c U(V)h u 
d v = i v G — d u (V)U — U(V)d u ~ ★ f — AU AV 



CO 

h v ( , TI ~ ~\ h u (V)~ U(V) 

— = i v *G = *(d u AUAV) + — —U 

c V / c c 
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Inserting these in the constitutive relations (|70l) yields 

d u (V)U - U(V)d u -*(^At/Ay)= £ [ e u (V)U - U(V)e u - * (c b u A U A Vj] 

c * (d u A U A + h Lr (V r )[/ - fj^h^ = ^ 1 [*(^A[/AV ; )+ b c/ (V r )[7 - CT(V)b u ] , 

where e = e r eo,fi = /i r /Xo in terms of the relative permittivity scalar e r and 
relative permeability [i r . Contracting with U gives the relations 



d u (V) = ee u (V) 
h u (V)= f i- 1 b u (V), 



(71) 



which yields 



U(V)d u + * ( — AUAV) = e t/(V> u + * ( cb AU AV 



U{V)b L 



• *i — At/Ay 

Co 



f7(F)h u - c * (d 17 A U A V^j 



(72) 



A laboratory Minkowski inertial frame is described by U — -^-dt in inertial 
Cartesian coordinates {t, x 1 , x 2 , x 3 } with {e° = —c dt, e 1 — dx 1 , e 2 = dx 2 ,e 3 = 
dx 3 } and, relative to [/, the medium 4- velocity X 7 has the orthogonal decompo- 
sition: 



V 



where the spatial field V' 7 is the Newtonian velocity fielcH of the medium relative 



tO U, EE 

that: 



and v 2 ee g(\Y u ,W u ). With these definitions it follows 



g(U,V) = -j 



and 



U A V = —UAY U . 



Co 



Using these in (|72[) and rearranging yields 

h^ 7 



d* 7 - * I 17 A Y u A 



u 



b + * [ £/ A V A — 

c o 



e^ 7 - * ( [7 A Y u A b u 



h u + * ( U A V u A d u 



(73) 



or using the spatial Hodge map # defined by U : 



d" + #(Y U A^\ =e e u +#(V c/ Ab c/ 



b — # V A 



(74) 



h u - # ( V 6 ' A d u 



These are the constitutive relations first written by Minkowski. For some pur- 
poses it is useful to decouple these expressions and express e u and h u directly 
in terms of d u , b u and V u . 



5 In inertial Cartesian coordinates V u = V' 3 , v J 'ft, x 1 , x 2 , x 3 ) -=2- . 
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D.2 Minkowski Constitutive Relations for e u (d u ,b u ,Y u ),h u (d u ,b u ,V U ) 
Taking the exterior product of ([71)1 with V u yields, 

fiW u Ah 11 = Y u A h u — V u A # (f u A ^ + /iV^ 7 A # (v y A d 17 ) . 

For any spatial 1-form a u with respect to U one has: 

V 17 A # (V 7 A a 17 ) = -¥ u A*(UAY u Aa u ^) =-*i v u (uAY u Aa u ^j 
= v 2 *(uAa u )~a u (W u )*(uAW u ) 
= -v 2 #a u +a u (V u )j0 u . (75) 

Thus 



1 



= W u Ah u + U -#e u - nu 2 #d u - -^(1 - en<$)d u (y u )#V u , 
c ec 

using the identity d c/ (V' 7 ) = ee u (Y u ), obtained by contracting ([7D with V u . 
Since ## = 1 



V / V / c ec 

Substituting this into the first relation of (|74p yields 



l-^jd^-(e M c2-l)f#(v^Ab 



Since 



JL = J_ 

e M TV 2 ' 

where J\f 2 = e r fi r is the square of the medium refractive index scalar, this may 
be written 

Furthermore from (1741) 



e# (V 17 A e 17 ) = # (V 17 A d u ) + (¥^p-W u - ^h c/ ) - e (b u (Y u )V u - i?h u ) 
= # (v u A d u ) - + ^h v - \ {e(M% - 1) b u (Y u )Y u , 

using h u {Y u ) = fih u (Y u ), obtained by contracting ^ with Y u . Substituting 
this into the second relation of (|74p yields 
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Thus, the constitutive relations can also be written 



2 



In the non-relativistic limit (to first order in — ) these constitutive relations 
become 



e u « — - ( 1 - — : : /-b' 

e e r V 6 

V — ]#(V t, Ad c ' 



(76) 



E Polarization and Magnetization 

The polarization 2-form II in spacetime is defined by 

II = G-e a F. (77) 
The second macroscopic Maxwell equation may then be written 
e Q d*F = j-d*U = j+j p , 

where 

j p = -d*U (78) 

will be called the electric polarization current 3-form. With respect to any 
observer frame U its orthogonal decomposition is 

IT = p u AU-*(— All] = p u AU--M u , (79) 

where M. u = #m u and we call 

p u — ii/II and =iu*H 

Co 

the spatial polarization 1-form and magnetization 1-form respectively relative 
to U. The Hodge dual of II has the decomposition 

m u ~ 1 

*n = *(p u AU) + A U = P u + — m u A U, (80) 

c c 

where P u = #p l/ . From (JSQ), ® , (|Z7J and O it follows 

d u = e e u +p u and h=n Q 1 h u + m u . (81) 

From ([78]). ((79)) and (JSTJ one finds 

j p = - d P u ~—AU = -dP u + UA£uI> u -^AU 
c c 

= -dP u + — (c CuP u - dm u ) A U 

Co 

= -dP u + - (P u ~ dm u ) AU. 
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Writing the orthogonal decomposition of j p with respect to U as 

Co 

it follows that 

J u 1 / • \ 

— = -iu3 P = — [P U ~ dm u ) and p u = -(ij, * Jp ) * [/ - -dP u . 
Co Co v / 

In the frame U, J p and p p denote the induced electric polarization current 
density spatial 2-form and induced polarization charge density spatial 3-form 
respectively. In a similar manner 

dU = dp u AU — —dM. u = dp u AU — —dM u + —U A £{jM. u 

c c c 

= — (c dp u + £uM u ) A U - — dM u . 

co c 

with the orthogonal decomposition j rn = dll = A U + p„ where 



-Wm = — c dp u + — M u and = * j m ) *U = dM L 

Co Co V c / c 

denote the induced magnetization charge current density spatial 2-form and in- 
duced magnetization charge density spatial 3-form respectively in terms of 
and M u . 



Bibliography 

[1] T. Dereli, J. Gratus and R. W. Tucker, The Covariant Description of 
Electromagnetically Polarizable Media, Phys. Lett. A, 361 (2006) 190-193 

[2] T. Dereli, J. Gratus and R.W. Tucker, New Perspectives on the Rel- 
evance of Gravitation for the Covariant Description of Electromagnetically 
Polarizable Media, J. Phys A: Math. Theor., 10 (2007) 5695-5715 

[3] M. A. Player, On the Dragging of the Plane of Polarization of Light Prop- 
agating in a Rotating Medium, Proc. Roy. Soc. A, 349 (1976) 441-445 

[4] J. Gotte, S M. Barnett and M. Padgett, On the Dragging of Light by 
a Rotating Medium, Proc. Roy. Soc. A, 463 (2007) 2185-2194 

[5] R. W. Tucker and T. J. Walton, An Intrinsic Approach to Forces in 
Magnetoelectric Media, II Nuovo Cimento C, 32 (2009) 205-229 



29 



